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Abstract 

We construct self-gravitating axially symmetric sphaleron solutions of the 3-1-1 dimen¬ 
sional Skyrme model coupled to Einstein gravity. The solutions are static and asymptotically 
flat, they are characterized by two integers n and m, where ±n are the winding numbers of 
the constituents and the second integer m defines type of the solution. These conhguration 
correspond to the chains of charge n Skyrmions and charge —n anti-Skyrmions placed along 
the axis of symmetry in alternating order. We investigate the dependency of the masses of the 
gravitating sphalerons on the gravitational coupling. We hnd new chains of self-gravitating 
|n| = 1 Skyrmions-anti-Skyrmions which emerge at some critical non-zero value of the gravi¬ 
tational coupling and do not have flat space limit. In contrast, the branches of self-gravitating 
|n| > 2 Skyrmion-anti-Skyrmion chains emerge from the corresponding flat space conhgura- 
tions. In both cases these branches merge at some maximal value of the effective gravitational 
coupling the branches of different type. The branch of gravitating Skyrmion-anti-Skyrmion 
pair extends all the way back to the limit of vanishing coupling constant where solutions 
approach the corresponding generalised Bartnik-McKinnon solutions. The upper branch of 
gravitating Skyrmion-anti-Skyrmion-Skyrmion chain exist up to some critical value of the 
gravitational coupling at which the chain becomes broken. We further find that for small 
values of the coupling constant on the upper branches, the solutions correspond to composite 
systems, consisting of a scaled inner Einstein-Yang-Mills solution and outer Skyrmions which 
are separating from the inner conhguration. 
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1 Introduction 


There are many nonlinear classical field theories in flat spacetime which admit topologically stable 
soliton solutions. These are particle-like, globally regular localised field configurations with finite 
energy. Interesting examples in d = 3-1-1 dimensions are monopoles in the Yang-Mills-Higgs model 
mm, the solitons of the Skyrme model [3] and knotted configurations of the Faddeev-Skyrme model 

[4] . The Skyrme model is nonlinear scalar 0(4) sigma-model, under certain assumptions it can 
be derived from the expansion of the QCD low energy effective Lagrangian in the large Nc limit 

[5] , then the topological charge of the multisoliton configuration is set into correspondence to the 
physical baryon number. Recently, some modifications of Skyrme model was proposed to approach 
the topological bound [5] preserving the topological properties of the corresponding solitons. 

Interestingly, in the flat space both the Yang-Mills-Higgs theory and the Skyrme model admit 
classical sphaleron solutions which represent monopole-anti-monopole pair in a static equilibrium 
[7], and Skyrmion-anti-Skyrmion solution [5], respectively. Further, there are axially symmetric 
generalizations of these solutions which represent chains of interpolating Skyrmion-anti-Skyrmions 
[S] and monopole-antimonopole chains [TO], each carrying topological charge n in alternating order. 
For an even number of constituents, these chains reside in the topologically trivial sectors of both 
models. When the number of constituents is odd, the configurations represent deformations of the 
axially symmetric multisoliton of degree n. 

When gravity is coupled to the matter field, this has significant effect on these solutions. It 
turns out there are hairy black hole solutions of the Einstein-Skyrme theory [II1II11I3]. His¬ 
torically, that was the first example of constructions of hairy black holes. These solutions are 
stable, asymptotically flat and possess a regular horizon, furthermore, they may be viewed as 
bound states of Skyrmions and Schwarzschild black holes [14] . Axially-symmetric static solutions 
of the Einstein-Skyrme model with topological charge two were studied in [15] . Recently, in |16] . 
self-gravitating BPS Skyrmions were applied in description of bulk properties of neutron stars. 

The globally regular gravitating Skyrmions in asymptotically fiat space were studied in [HIS]. 
It was shown that there are two branches of solutions, one of which emerges smoothly from the fiat 
space Skyrmion configuration. When the effective gravitational coupling constant = dTrGF^, 
where G is the Newton constant and Fn- is the pion decay constant, is increased from zero, this 
branch terminates at some critical value of the coupling, beyond which gravity becomes too strong 
for self-gravitating Skyrmions to persist. There it merges with a second branch, which extends 
all the way back to vanishing coupling constant. Along this branch the mass of the gravitating 
Skyrmion rapidly increases and the solution becomes unstable. Surprisingly, it was shown [13] 
that in the limit of vanishing coupling, the gravitating Skyrmion approaches the lowest Bartnik- 
McKinnon (BM) solution of the SU{2) Einstein-Yang-Mills theory [17]. This pattern is rather 
similar to the branch structure of the gravitating monopole-antimonopole pair [18] . Indeed, the 
effective coupling constant of the Einstein-Yang-Mills-Higgs theory is proportional to the Higgs 
vacuum expectation value r] and the square root of the gravitational constant G. The evolution of 
the sphaleron solution along the second branch may be considered as being obtained by decreasing 
the vacuum expectation value of the scalar field, in the limit of 77 —>■ 0 the configuration smoothly 
approaches the same lowest mass spherically symmetric BM solution. 

The pattern of evolution of the self-gravitating monopole-anti-monopole chains and vortex rings 
in Einstein-Yang-Mills-Higgs theory is qualitatively similar to that of the monopole-antimonopole 
pair, the difference is that on the upper branches, these solutions correspond to composite systems, 
consisting of a scaled inner Einstein-Yang-Mills solution and an outer flat space Yang-Mills-Higgs 
solution m- 

The properties of gravitating Skyrmions were considered in many works, for example config¬ 
urations with discrete symmetry were investigated in |20j . spinning gravitating Skyrmions were 
studied in |2T]. Also modifications of the Einstein-Skyrme model with cosmological constant were 
investigated in |221123] . In this paper we present globally regular gravitating axially-symmetric 
Skyrmion-anti-Skyrmion configurations. 

Since the consistent consideration of the solitons with higher number of constituents is related 
with complicated task of numerical simulations, we restrict our consideration to the case of the 


2 


Skyrmion-anti-Skyrmion (S-A) pair and Skyrmion-anti-Skyrmion-Skyrmion (S-A-S) chain with 
constituents of degrees n = 1,2. We show that the general pattern of evolution of the configuration 
is similar to the branch structure of the monopole-antimonopole chains in the Einstein-Yang-Mills- 
Higgs system which links the corresponding flat space configurations and the BM solution. 


2 The model 


The Einstein-Skyrme model in asymptotically flat 3-1-1 dimensional space is defined by the action 


S = 


R 


IGttG 


+ Lm > 


( 1 ) 


where the gravity part of the action is the usual Einstein-Hilbert action with curvature scalar R, 
g denotes the determinant of the metric, G is the gravitational constant and the matter part of 
the action is given by the Skyrme Lagrangian 


= ^gf^-Tv{L^L, 




La-]) -L 




Tr {U - 1) 


( 2 ) 


where = V^d^U is the su(2)-valued left current, associated with the S'{7(2)-valued Skyrme 
field U = cr • I -I- i 7 r“ • t“. Here the quartet of the fields (o’, 7 r“) is restricted to the surface of the 
unit sphere, cr^ -|- 7 r“ • 7 r“ = 1 and e and are parameters of the Skyrme model. Explicitly, 
En- is the so-called pion decay constant, e is dimensionless constant and is the tree-level pion 
mass. Note that in the rescaled radial coordinate r —>■ erEn./2 the action of the Einstein-Skyrme 
theory becomes 


^ = j |^ + ^<?'^"Tr(E^E.) + l5'^V"Tr([E^,E,][E,.,E<,])+M"Tr(E-l)|^d4 


(3) 


where g = 2^,.,^/[FT^e) is the rescaled mass parameter. Thus, the only physical parameter of the 
Einstein-Skyrme theory ([3|) is the effective gravitational coupling constant = AttGFt^ [T3]. 
Variation of the action m with respect to the metric 5 ^ 1 , yields the rescaled Einstein equations 




(4) 


where Skyrme stress-energy tensor is 


i;. =Tr - L^L^j + TrE^][E“, 

- 45„^[E^, E“][E„ E^]^ + M' 5 M.Tr(t/ - 1). 


(5) 


In order to construct finite energy topologically non-trivial solutions of the model ([3]) the 
Skyrme field must approach the unit matrix at spacial boundary, E —>■ 1 as r —>■ 00 , thus the field 
is a map U : ^ S^, which belongs to an equivalence class characterized by the homotopy group 

tt3{SU{2)) = Z. The corresponding topological current is 

= -_^=e^,^,Tr(E'^E^E'^). (6) 

To obtain gravitating static axially symmetric solutions, we employ the usual Lewis-Papapetrou 
metric in isotropic coordinates: 


ds^ = -fdi‘‘ + + j 


dO^ 


Ir^ 


sm 


-d(p'^ 


(7) 
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where the metric functions /, m and I are functions of the radial variable r and polar angle 0, 
only. The z-axis (0 = 0,7r) represents the symmetry axis. 

Since we are only concerned with axially symmetric fields, the SU{2) valued Skyrme field 
U{r) = cr • I + i7r“ • r“ can be parameterized as [SI H] 

TT^ = (j)^ cos{n(p); = (j)^ sm{rup)-, tt^ = a = (j)^ (8) 

where the triplet of scalar fields </)“ on unit sphere is a function only of radial variable r and polar 
angle 9. One can check that the parametrization ([8]) is consistent, i.e. the complete set of the 
field equations, which follows from the variation of the action of the Einstein-Skyrme model, is 
compatible with the equations on three functions (/>“, which follow from variation of the reduced 
action on the ansatz (jS)). 

An equivalent parametrization of the axially-symmetric Skyrme field is given by 

(j)^ = sin f sing] = sin/cosg; (j)^=cosf (9) 

where the functions f{r,9) and g{r,9) depend on the radial coordinate r and polar angle 9. 
Imposing the boundary condition on the asymptotic value of the function g{r,9) [S] 

lim g{r,9) = m9 (10) 

r—^oo 

where m is an integer, we can construct various multi-Skyrmion or the Skyrmion-anti-Skyrmion 
chain solution of given degree. Then the value of m corresponds to the number of the constituents 
of the configuration, each of them can be identified with an individual charge ±n Skyrmion. 
Evaluation of the corresponding baryon number yields B = f [1 — (—1)™], i.e. the case m = 
1 corresponds to the multi-Skyrmions, in particular n = 1 ,to = 1 configuration is the usual 
spherically symmetric Skyrmion parameterized via the function of radial variable f{r). 

The m = 2 configuration is the Skyrmion-anti-Skyrmion (S-A) pair |8], which is a deformation 
of the topologically trivial sector. Note that in the absence of the potential term in the flat space 
only topological charge pairs of n > 2 are found. As we will see, the gravitational interaction 
between the constituents makes it possible for n = 1 self-gravitating Skyrmion-anti-Skyrmion 
pairs to exist. 

The configurations with m = 3 correspond to the Skyrmion-anti-Skyrmion-Skyrmion chain (S- 
A-S) M- Thereafter we restrict our consideration to the axially symmetric solutions with values 
n = 1, 2 and m = 1, 2, 3. Note that the binding energy of the constituents of the Skyrmion-anti- 
Skyrmion chains is quite weak, so the \n\ = 2 chains in the flat space are very unstable with 
respect to perturbations [9]. 

Clearly, these configurations and monopole-anti-monopole chain solutions of the Yang-Mills- 
Higgs theory laiio] have many features in common. In both cases there is an effective interaction 
between the constituents which allows the sphaleron solution to exist, although the nature of the 
interaction is different. In the Yang-Mills-Higgs theory it is an effective electromagnetic interaction 
between the monopoles [24] while in the case of the Skyrme theory there is a dipole-dipole inter¬ 
action between the Skyrmions [SlIS]. Since the latter interaction is much weaker than the effective 
electromagnetic forces, in the absence of the mass term the Skyrmion-anti-Skyrmion chains may 
exist in the flat space only when each of the constituents carries charge |n| > 2. Clearly, cou¬ 
pling to gravity provides additional attraction between the lumps which results in the existence of 
self-gravitating Skyrmion-anti-Skyrmion chains with |n| = 1. Thereafter we set the mass constant 

fj, = 0. 

The complete set of the field equations, which can be obtained from variation of the action of 
the Einstein-Skyrme model (jS)), can be solved when we impose the boundary conditions on the 
fields (jf" and make use of the parametrization of the metric 0. Then the field equations reduce 
to a set of six coupled partial differential equations, to be solved numerically. 

As usually, the boundary conditions follow from the regularity on the symmetry axis and 
symmetry requirements as well as the condition of finiteness of the energy. In particular we have 
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to take into account that the asymptotic value of the Skyrme field is restricted to the vacuum and 
the metric functions must approach unity at the spacial boundary. Explicitly, we impose 


0 , 

1, m 


■ 0 , 

1 , I 




at infinity. For odd values of m at the origin we impose 


0 , 


drf 


r^O 


r —>-0 


0 , 


1 —^0 


- 1 . 


r —>-0 


0 , drUl 


0 , dJ 


r —)-0 


^0. 


r —>-0 


For the S-A pair the boundary conditions on the Skyrme field are different 019] 


( 11 ) 


( 12 ) 


0 , 


-->-0 


0 , 


r —>-0 


r —)-0 


(13) 


since the position of the components corresponds to the points in space where the third component 
of the field (j)^ is equal to its anti-vacuum value —1. 
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Figure 1: The mass M (left), and the scaled mass M (right) of the charge 1 gravitating Skyrmion and the 
Skyrmion-anti-Skyrmion chains for n = 1 , m = 1, 2, 3 are shown as functions of the coupling constant a. 
The thin lines extend the Skyrmion curves of the scaled mass to the mass of the corresponding generalized 
Bartnik-McKinnon solution. 


The condition of regularity of the functions on the symmetry axis yields 



—)■ 0, 


—)■ 0, 

dgc, 

(>3 

—^ 0 



9 ^ 0,77 


>-0,7r 



9—^0,7t 

(14) 

def 

—)■ 0 , 

dgm 

—)■ 0, 

dgl 


0 


0—>-0,7r 


0—f0,7r 


9 - 

-4-0,77 



To satisfy the condition of regularity m(r, 0) = Z(r, 0) we introduce the auxiliary function g{r, 0) = 
l{r,9)/m{r,9) with the boundary conditions g{0,9) = g{oo,9) — g{r,0) = 1; dgg{r,Trl2) = 0). We 
check this condition as a test for correctness of our numerical results. 


3 Numerical results 

The numerical calculations are mainly performed on an equidistant grid in spherical coordinates 
r and 9, employing the compact radial coordinate x = r/(l -I- r) G [0:1] and 9 G [0,7r]. To find 
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solutions of the Euler-Lagrange equations which follow from the rescaled action (jS]) and depend 
parametrically on the effective gravity constant a, we used the software package CADSOL based 
on the Newton-Raphson algorithm [25) . This code solves a given system of nonlinear partial 
differential equations subject to a set of boundary conditions on a rectangular domain. Typical 
grids we used have sizes 90 x 70. The resulting system is solved iteratively until convergence is 
achieved. 

Apart from some initial guess for the solution, CADSOL requires also the Jacobian matrices 
for the equations with respect to the unknown functions and their first and second derivatives, and 
the boundary conditions. This software package provides also error estimates for each function, 
which allows to judge the quality of the computed solution. The relative errors of the solutions we 
found, are of order of 10“^ except the S-A-S chains for those they are as large as 10“^. We also 
introduce an additional Lagrangian multiplier to constrain the field to the surface of unit sphere. 
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Figure 2: The mass M (left), and the scaled mass M (right) of the charge 2 gravitating Skyrmion and the 
Skyrmion-anti-Skyrmion chains for n = 2 , m = 1, 2, 3 are shown as functions of the coupling constant a. 
The thin lines extend the Skyrmion curves of the scaled mass to the mass of the corresponding generalized 
Bartnik-McKinnon solution. 




Figure 3: The values of the metric functions /(O) at the origin are shown as functions of the gravitational 
coupling constant a for the gravitating Skyrmion and the Skyrmion-anti-Skyrmion chains with n = 1 
(left) and n = 2 (right) with m = 1, 2, 3. 

Note that dimensionless gravitational coupling constant = AttGFt^ vanishes if (i) the Newton 
constant G —>■ 0, or, (ii) JG —>■ 0. In the former case one may recover the sphaleron solitons of 
the Skyrme model in the flat space, in the second case the sigma-model term in the action o 
is vanishing. Thus, similar to the case of the self-gravitating monopole-antimonopole systems in 
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the asymptotically flat Einstein-Yang-Mills-Higgs theory mns] and Yang-Mills sphalerons in the 
AdS 4 spacetime [26) . there are two branches of solutions of the Einstein-Skyrme model [HHZlIlI]. 

The dependence of the spherically symmetric n = l,m = 1 Skyrmion on gravity has been 
studied before [nmnii, this branch of gravitating Skyrmions is linked to the flat space solu¬ 
tions. Similarly, when gravity is coupled to the Skyrme model, a branch of gravitating axially- 
symmetric Skyrmion-anti-Skyrmion chains, which are formed from m constituents of charge 2, 
emerges smoothly from the corresponding flat space sphaleron solutions. 

However, |n| = 1 sphalerons of the massless Skyrme model do not exist in the flat space [ 8 ]. 
Clearly, coupling to gravity yields additional attractive interaction between the constituents. We 
have found numerical evidence that in the Einstein-Skyrme model a branch of self-gravitating 
charge 1 Skyrmion and charge -1 anti-Skyrmion pairs emerges at the critical non-zero value of 
the gravitational coupling ao = 0.08 and extends up to a maximal value ttcr = 0.251 where it 
merges the upper mass branch. For the \n\ = 1 S-A-S chain the critical values of the gravitational 
coupling are ao = 0.082 and acr = 0.2305, respectively. Along this so-called Skyrmion branch, the 
distance between the constituents, which is defined as the points in space where the component 
of the Skyrme field (^3 is equal to its anti-vacuum value —1, is slowly decreasing. 

As the gravitational coupling constant increases, the background is getting more and more 
deformed and, at some critical value of the coupling, when the gravity becomes too strong for 
solutions to persist, the Skyrmion branch merge with branch of different type, cf Figs. [l][31 The 
critical value acr, at which a backbending is observed, slightly decreases as the number of the 
constituents of the Skyrmion-anti-Skyrmion chains is increasing. 

The parametrization m allows us to find the dimensionless ADM mass of the configuration 
M, it is defined by the value of the derivative of the metric function / at the boundary 

M = lim da;f (15) 

x—¥oo 

In order to perform another check of our numerics for correctness, we compare this value with the 
results of direct evaluation of the integral over the Tqq component of the total energy-momentum 
tensor ®. 

Along the first (lower) branch the mass of the gravitating Skyrmion-anti-Skyrmion chains 
decreases with increasing a, since with increasing of gravitational strength the attraction in the 
system increases. Along the second (upper) branch, in contrast, the mass p5p increases strongly 
with decreasing the coupling a, and the solutions shrink correspondingly. However, in the limit 
of vanishing coupling constant when the mass M diverges, the pattern of evolution of the config¬ 
urations depends on the value of integer numbers n and m. 

As shown in m the spherically symmetric n = 1 , to = 1 configuration in this limit approaches 
the lowest mass spherically symmetric Bartnik-McKinnon solution m- The structure of the 
limiting ^ 0 configuration can be better understood when we introduce the rescaled radial 
coordinate x = xj a and the rescaled mass M = aM [T3]. Then the reparametrization of the 
Skyrmion field 

4>^{x) = \/l — uj[xY sin 0 ; = ■\/l — oj{xY cos 0 ; 4>^{x)=uj{x) (16) 

allows us to represent the system of the field equations on the function u)(x) and the metric 
functions f{x),l{x) = m{x) in the form, which is identical to the system of the SU{2) Einstein- 
Yang-Mills equations for the Bartnik-McKinnon configuration with the gauge connection A“ = 
eian(l — oj{x))Xn/2 [TT] . On the other hand, the rescaled mass of the axially symmetric charge 
n = 2 Skyrmion is higher, the limiting configuration on the corresponding upper branch is the 
generalized BM solution [28l [29] . 

Let us now consider axially-symmetric to = 2 Skyrmion-anti-Skyrmion pairs. This configura¬ 
tion resides in the vacuum sector, these chains exhibit an analogous dependence on the coupling 
constant a as the to = 1 Skyrmions. We found that this self-gravitating n = 1 and n = 2 S-A 
sphaleron with decreasing a on the upper branch also smoothly approaches the corresponding 
generalized BM solution, as illustrated in Figs. I1I2I 
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To clarify this observation, we plot in FigUjthe metric function / and the third component of 
the Skyrme field for the n = 1 conhgurations for a small value of the gravitational coupling a on 
the upper branches. 

Clearly, we can identify three distinct regions. As seen in Fig. 31 right column, in the first 
region the metric function / remains very small but constant without any significant angular 
dependency. In the second transition region, the metric varies up to upper value / = 1, this is a 
small region where the energy of the matter field of the S-A configuration is located. Finally, in 
the third outer region, the metric functions are approaching the flat space limiting values. 

This pattern is similar to the known picture of the evolution of self-gravitating monopole- 
antiminopole pairs mus], which on the upper unstable branch are also linked to the generalized 
Bartnik-McKinnon limiting solutions. One might expect the same analogy holds for the Skyrmion- 
anti-Skyrmion gravitating chains and monopole-antimonopole chains with higher number of con¬ 
stituents. However, we observe a bit different pattern for the coupling constant dependence of the 
S-A-S system on the upper branch. For some small critical value of the gravitational constant a 
the central component of the chain, the anti-Skyrmion of charge —n, is located into the interior 
region while remaining in the transition region pair of Skyrmions with positive charges n becomes 
separated since the weak gravitational interaction there cannot stabilize it. Thus, the n = 1,2 
S-A-S chains are broken at = 0.06 and = 0.075, respectively. 

The composition of the n = 2 ,to = 3 solution at — 0.01 is exhibited in Fig. [SJ in Fig. 31 
bottom row, we presented the the metric function / and the third component of the Skyrme held 
for the n = 1 S-A-S conhguration. 

However, one may expect that for higher values of the components of the chains, say for m = A 
S-A-S-A conhguration with \n\ > 3, will stay unbroken along the upper branch all the way down 
to the limit a —>■ 0. This conjecture is based on the preliminary results we found considering 
TO = 4, n = 3 solution, they suggest that as a decreases, the location of the inner S-A pair moves 
continuously inwards along this branch and tends to the origin in the limit a —>■ 0. However, 
the attractive interaction between the remaining in the outer hat space region S-A pair could be 
strong enough to overcome the repulsive force. The location of the outer S-A pair in this case 
approaches a hnite value, which is in a good agreement with the separation of the corresponding 
S-A pair in the hat space. Thus, similar to the case of to = 4 monopole-antimonopole chain, the 
resulting conhguration may be thought of as composed of a scaled (generalized) BM solution in 
the inner region and a hat space soliton-anti-soliton pair solution in the outer region. 


4 Conclusions 

The main purpose of this work was to present new type of self-gravitating sphaleron solitons in 
the Einstein-Skyrme theory. The conhgurations of that type are characterized by two integers n 
and TO, where \n\ is the winding numbers of the constituents, and the second integer to dehnes 
the type of the conhguration, it is a deformation of the topologically trivial secror for even to and 
for odd values of to the Skyrmionanti-Skyrmion chain is a deformation of the axially-symmetric 
Skyrmion of degree n [9]. 

In the present work we have focussed on the Skyrmion-anti-Skyrmion chains with 2 and 3 
constituents, S-A and S-A-S, respectively. Winding number of each individual component is 
restricted to lower values n = 1,2. These solutions are asymptotically hat, and globally regular. 

Concerning the dependence of the gravitating Skyrmion-anti-Skyrmion chains on the gravity 
coupling constant, we generally observe the pattern which is similar yet slightly different from the 
picture observed previously in the case of the sphaleron solutions of the Einstein-Yang-Mills-Higgs 
theory. Eirst, in the absence of the pion-mass potential, there are branches of self-gravitating \n\ = 
1 Skyrmion-anti-Skyrmion chains which emerge at some critical non-zero value of the gravitational 
coupling, these branches do not have hat space limit. In contrast, the branches of self-gravitating 
\n\ > 2 Skyrmion-anti-Skyrmion chains emerge from the corresponding hat space conhgurations. 

In both case these lower branches merge at some maximal value of the effective gravitational 
coupling the upper branches. The pattern of evolution along these branches is related with effective 
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decreasing of the Skyrme coupling constant and the structure of the limiting configurations 
depends on the type of the configuration. The branches of gravitating m = 1 Skyrmions and m = 2 
S-A pairs extends all the way back to the limit a —> 0 where solutions approach the corresponding 
(generalised) Bartnik-McKinnon solutions of the SU{2) Einstein-Yang-Mills theory. However, the 
upper branch of gravitating S-A-S chain exist up to some critical value of the gravitational coupling 
at which the chain becomes broken. Then the anti-Skyrmion of degree n remains constrained at 
the interior region whereas the pair of identically charged Skyrmions in the outer region becomes 
separated. 

Here it is interesting to compare the results with the pattern of evolution of the self-gravitating 
solitons of the Einstein-Faddeev-Skyrme model [30]. The structure of the Skyrme model and the 
Faddeev-Skyrme model look similar, the corresponding Lagrangians, like ([2), include the usual 
sigma model term, the Skyrme term, which is quartic in derivatives of the scalar held (jf", and 
the optional potential term which does not contain the derivatives. However the topological 
properties of the corresponding solitons are different, the hnite energy solutions of the Faddeev- 
Skyrme model, the HopHons, correspond to the map <()“ : —?> 5^ which belongs to an equivalence 
class characterized by the third homotopy group = Z [4]. Thus, the corresponding integer 

topological invariant associated with the triplet of scalar helds (()“ constrained to the unit sphere 
S'^ is known as the Hopf number Q. This invariant can be interpreted geometrically as the linking 
number of two loops obtained as the preimages of any two generic distinct points on the target 
space S^. 

Similar to the case of the Einstein-Skyrme model considered above, we can couple the scalar 
triplet of the fields of the Faddeev-Skyrme model to gravity. Considering the gravitating static 
axially symmetric Hopfions of lower degree Q < 4 we observe the same general pattern as for to = 1 
Skyrmion solutions of the Einstein-Skyrme model, S-A chains and the monopole-antimonopole 
chains of the Einstein-Yang-Mills-Higgs theory. We have found numerical evidence that, when 
gravity is coupled to the Fadeev-Skyrme model, a branch of gravitating Hopfions emerges from 
the flat space Hopfion solution and extends up to a maximal value acr where it merges the upper 
mass branch, cf Fig. |6| vs Figs. I1I2I Once again, the upper branch extends back to the limit 
a —^ 0 where solutions approach the corresponding (generalised) Bartnik-McKinnon solutions of 
the SU{2) Einstein-Yang-Mills theory, as seen in Fig..|SJ right plot. 

We hope to return elsewhere with a discussion of some of these interesting aspects of the 
self-gravitating Hopfions. 

There are various possible extensions of the solutions discussed in this work. Clearly, it would 
be interesting to investigate how inclusion of the pion mass term may affect the results. Also, we 
expect the to > 4 Skirmion-anti-Skyrmion chain configurations with constituents of higher degree 
\n\ > 3 may stay unbroken on the upper branch as a —> 0. Finally, our preliminary results indicate 
the existence of static axially symmetric black hole solutions with Skyrmion’s hair, consisting of 
chain-like structures. It would be also interesting to address the question how inclusions of a 
cosmological constant will affect the properties of a gravitating sphalerons of the Skyrme model. 
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Figure 4: The metric function /(r, 6) (right column) and the third components of the field of the gravitat¬ 
ing spherically symmetric Skyrmion and the Skyrmion-anti-Skyrmion chains with components of charge 
n = 1 (left column) on the upper branch are shown in logarithmic radial coordinate logr at o? — 0.009 
(for m = 1, 2) and at o? — 0.01 (for m = 3. 


11 













































Figure 5: The components of the field of the gravitating n = 2,m = 3 S-A-S chain on the upper branch 
are shown at c? — 0.01 in the logarithmic radial coordinate logr. 




Figure 6: The mass and the scaled mass fi of the gravitating Hopfions of degrees 1 to 4 are shown as 
functions of the coupling constant a. The thin lines extend the Hopfion curves of the scaled mass to the 
mass of the corresponding generalized Bartnik-McKinnon solution. 


12 


















